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Abstract. We present deep V and I photometry for the 
open cluster NGC 4815 and four surrounding Galactic 
fields down to a limiting magnitude V ~ 25. These data 
are used to study cluster spatial extension by means of star 
counts, and to derive the luminosity (LF) and mass func- 
tion (MF). The radius turns out to be 3.6 ± 0.3arcrran at 
V=19.0 level, whereas the mass amounts at 880 ± 23Om,0 
down to V=20.8. 

From the color-magnitude diagram, we obtain the LFs in 
the V and I bands, using both the standard histogram and 
an adaptive kernel. After correction for incompleteness 
and field star contamination, the LFs have been trans- 
formed into the present day mass functions (PDMF). The 
PDMFs from the V and I photometry can be represented 
as a power-law with a slope a = 3.1±0.3 and a — 2.9±0.3 
(the (?) MF in this notation has a slope a = 2.35) respec- 
tively, in the mass range 2.5 < ^L. < 0.8. 
Below this mass, the MF cannot be considered as rep- 
resentative of the cluster IMF, as it is the result of the 
combined effect of strong irregularities in the stellar back- 
ground, probable internal dynamical evolution of the clus- 
ter and/or interaction of the cluster with the dense Galac- 
tic field. Unresolved binaries and mass segregation can 
only flatten the apparent derived IMF, so we expect that 
the real IMF must be steeper than the quoted slope by an 
unknown amount. 
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1. Introduction 



Galactic open clusters are loose and weakly bound 
systems, which only in exceptional cases survive more 
than 10 8 yrs (?). Their mean lifetime is determined 
essentially by two factors: the birth-place (or the orbit) 
and the initial mass. The dissolution of open star clusters 
is mainly due to close encounters with giant molecular 
clouds which are very frequent in the galactic thin disk 
(?; ?). For this reason, the age distribution is strongly 
biased towards young objects (?). Out of about 1,200 
open clusters which are known in the thin disk, only 80 
are older than the Hyades (?). 

The dynamical evolution of open star clusters in the 
Galactic disk environment has been studied by several 
authors (?; ?; ?) To summarize, it has been shown that 
encounters with dense molecular clouds are catastrophic, 
whereas diffuse clouds do not completely destroy the 
cluster. Moreover, tidal effects induced by the general 
Galactic gravitational field disturb significantly cluster 
evolution, producing distorted debris, at the center of 
which a few closely bound stars remain, with a large 
fraction of binaries. 

On the other hand, the internal dynamics of open clusters 
has been investigated by means of N-body models (?; ?), 
which include the effects of stellar evolution and a fine 
treatment of binaries and/or multiple stars. 
Among the others, a very important result is that the 
cluster evolution crucially depends on the Initial Mass 
Function (IMF), and on the cluster richness. To go in 
some detail, de la Fuente Marcos (?) finds that power-law 
(?) like IMFs accelerate the cluster disruption if the 
initial population is small, and that at the end of the 
evolution a huge number of binary systems with the same 
mass is formed. 
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An important constraint for all these models would 
be the possibility to derive the IMF from observations of 
open star clusters. Unfortunately, this is a very difficult 
task for several reasons. 

First, most open clusters are poorly populated ensembles 
of star, containing from tens to hundreds of stars, and 
only in rare cases up to some thousands. Secondly, they 
are located in the inner regions of the thin disk. This fact 
renders it very difficult to detect cluster members (at odds 
with globular clusters), due to the high contamination of 
the field stars in the same line of sight as the clusters. 
Third, they can be strongly obscured due to interstellar 
absorption lying between us and the cluster (?). 

Probably, the most challenging task is to obtain clus- 
ter membership. This is feasible for objects close to the 
Sun, for which we can obtain high quality radial veloc- 
ity and proper motion measurements. Going farther away 
one has to rely on statistical corrections which usually 
are made by comparing the cluster with the field in its 
outskirts. Ideally, when a correction for field star contam- 
ination has been performed, a reasonable actual (Present 
Day) mass function (PDMF) can be derived from the ob- 
served LF. By modeling the dynamical evolution of the 
cluster the PDMF can be converted into an IMF which 
can be compared with the MFs of other clusters to look 
for universality or deviations. 

In this paper, the first of a series dedicated to obtaining 
luminosity and mass functions for Galactic open clusters, 
we present results for NGC 4815, a rather rich, distant 
open cluster located at about 2.5 kpc from the Sun . 
The ultimate goal is to provide a sample of PDMFs as a 
function of cluster age, population, and location to serve 
as template for N-body models. 

We plan to study other clusters with the same method in 
forthcoming papers. In order to obtain a good representa- 
tion of the field stellar population, in addition to the clus- 
ter we observed four nearby Galactic fields, whose popu- 
lation has been averaged and used for the statistical sub- 
traction of the foreground/background objects from the 
cluster. An adaptive kernel estimator was then applied 
to obtain a continuous LF, free of binning bias, together 
with the standard histogram technique. Finally, the MF 
has been derived from m — AmriQ down to m — O.35m0 by 
using the luminosity-mass relation from the Girardi et al. 
(?) models. 

The layout of the paper is as follows. Section 2 summa- 
rizes our present knowledge of the open cluster NGC 4815. 
In Section 3 we present the observations and data reduc- 
tion technique; Section 4 is dedicated to the discussion of 
the CMDs obtained for the cluster and the field. Section 5 
is dedicated to a detailed analysis of the star counts, while 
Section 6 deals with the LFs derived with the two differ- 
ent techniques. Section 7 presents the MF of NGC 4815 
and the comparison with other MFs. Finally Section 8 
summarizes our results. 



Table 1. NGC 4815 fundamental parameters. Age, dis- 
tance and reddening are taken from the study of Carraro 
and Ortolani (1994) 



I 
b 

E(B - V) 

(to - M) v 

(m-M)j 

Age 

Z 



303.63 
-2.09 

0.70±0.05 

14.10±0.05 

13.20±0.05 

(0.5 ±0.25) x 10V 

«0.008 



2. Previous investigations 

NGC 4815 has been studied in the past by several authors. 
The first photometric data have been collected by Moffat 
& Vogt (?), who studied 9 stars in the cluster region, but 
did not recognize any sequence. 

Kjeldsen & Frandsen (?) presented a catalogue of UBV 
photometry for 559 stars. 

Carraro and Ortolani (?) firstly made an accurate analy- 
sis of the CMD by obtaining BV photometry of 2500 stars 
down to V=19.5. Their results led to the determination 
of the cluster parameters: NGC 4815 turns out to be of 
intermediate age (about half a Gyr old, like the Hyades) , 
with a reddening E(B — V) ~ 0.70, and a distance of 
about 2.5 kpc from the Sun. The metallicity is half Solar 
(see Table 1). From a preliminary analysis of the MS mass 
function they obtained an index a = 2.55, slightly higher 
than the classical Salpeter value. 

These data have been used to study the cluster struc- 
ture and luminosity function with a different method by 
Chen et al. (?), who obtained for the first time the clus- 
ter mass function by subtracting the field contamina- 
tion with a Galaxy structure model. In the mass range 
1.1 < — < 2.2 they obtain a slope a = 1.97±0.17, and a 
total mass of 912 m Q to the photometric limit (V = 19). 
They finally found evidence for mass segregation. 

3. Observations and reduction 

3.1. Observations 

NGC 4815 has been observed with the ESO-NTT tele- 
scope on April 14-16, 1999. We used the EMMI red arm, 
equipped with the 2048x2048 ESO#36 Tektronics CCD. 
The scale on the sky is 0.27 arcsec/pixel, for a total field 
of 9.2 x 9.2 arcmin 2 . The observing log is in Table 1. 
Almost all the images have been collected in sub-arcsec 
seeing conditions and in photometric conditions. During 
the same nights, four Galactic fields have been observed. 
These fields are located at 30 arcmin from the center of 
NGC 4815, on the North, South, East and West sides, re- 
spectively. The exposure times and the seeing of the Galac- 
tic field images are the same as for the cluster data, with 
the only exception of the southern one, for which only 
short exposures could be obtained. We opted to exclude 
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this field in the statistical analysis of luminosity and mass 
function. 

Table 2. Log-book of the NGC 4815 observations 

Field Filter Exp. Time) Seeing 

[s] FWHM["] 


15 




Center 


V 


15 


0.8 






Center 


I 


15 


0.8 






Center 


V 


900 


0.8 






Center 


I 


900 


0.8 


20 




West 


V 


15 


1.0 


> 




West 


I 


15 


1.1 






West 


V 


900 


1.0 






West 


I 


900 


0.9 






North 


V 


15 


0.7 






North 


I 


15 


0.8 


25 




North 


V 


900 


1.0 






North 


I 


900 


0.9 






East 


V 


15 


0.7 






East 


I 


15 


0.7 






East 


V 


900 


0.9 






East 


I 


900 


0.9 


30 


2 


South 


V 


15 


0.8 






South 


I 


15 


0.6 







3.2. Data reduction and calibration 

All the images have been pre-processed in the standard 
way with IRAF, using the sets of bias and sky flat- 
field images collected during the observing nights. The 
instrumental magnitudes, and the positions of the stars 
for each frame have been derived by profile-fitting pho- 
tometry with the package DAOPHOT II and ALLSTAR 
(?). Then we used ALLFRAME (?) to obtain the si- 
multaneous PSF-fitting photometry of all the individual 
frames. In order to obtain the transformation equations 
relating the instrumental (vi) magnitudes to the stan- 
dard V (Johnson), /(Kron-Cousins) system we followed 
the procedure already described in (?). Seven Landolt 
(1982) fields of standards have been observed. Specifically: 
the PG0918+029, PG0942-029, PG1047+003, PG1525- 
071, PG1657+078, SA104334, and SA110-496 regions. In 
each of these, there are other secondary standard stars by 
Stetson (?), which extend the previous Landolt sequence. 
For all of these stars, aperture photometry was obtained 
on all the images after the digital subtraction of neigh- 
boring objects from the frames. We used DAOGROW (?) 
to obtain the aperture growth curve for each frame and 
other related codes (by PBS) to determine the aperture 
corrections, the transformation coefficients, and the final 
calibrated photometry (e.g., Stetson 1993). We used trans- 
formation equations of the form: 

v = V + A + A x * X + A 2 * (V - I), 




Fig. 1. The CMD of NGC 4815. A selection based on the 
sharp and chi parameters has been done following Stetson 
(1987). Sharp is vaguely related to the intrinsic (i. e. out- 
side the atmosphere) angular size of the astronomical ob- 
ject. For an isolated star, sharp should have a value close to 
zero, whereas for semi-resolved galaxies and unrecognized 
blended doubles sharp will be significantly greater than 
zero, and for cosmic rays and some image defects will be 
significantly less than zero. Chi is a robust estimate of the 
observed pixel-to-pixel scatter from the model image pro- 
file divided by the expected pixel-to-pixel scatter from the 
image profile. The straight lines indicate the CMD region 
used to derive the LFs and represents the limiting magni- 
tude for the 50% completeness level. Finally, the crosses 
on the left side of the plot represent the photometric errors 
at that magnitude level. 

i = I + Bo + B x * X + B 2 * (V - /), 



Table 3. Coefficients of the calibration equations 





April 14,1999 


April 15,1999 


A 


0.053±0.007 


0.044±0.006 


Ai 


0.249±0.041 


0.259±0.023 


M 


-0.036±0.002 


-0.028±0.003 


B 


0.690±0.012 


0.601±0.012 


Si 


-0.054±0.061 


0.307±0.050 


B 2 


0.020±0.002 


0.030±0.005 



where the values of the coefficients are listed in Table 3 for 
the two nights. In these equations v and i are the aperture 
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magnitude already normalized to 1 sec, and X is the air- 
mass. Second order color terms were tried and turned out 
to be negligible in comparison to their uncertainties. It 
is a reasonable hypothesis that the color-dependent terms 
imposed by the telescope, filter, and detector should be 
reasonably constant over the course of a few nights, so af- 
ter A2 and B2 had been independently determined from 
each night's data, we obtained weighted mean values of 
these parameters and imposed them as fixed quantities 
for the definitive calibration. The NGC4815 images en- 
abled us to define local standards that met the follow- 
ing conditions: each was well separated from its neighbors 
(according to a criterion described in Stetson (?) §4.1), 
each was observed in all frames, each had a mean value of 
the goodness-of-fit statistic, CHI, less than 1.5. Once cal- 
ibrated, this local standard sequence enables us to place 
all of observations of NGC 4815, and of the Galactic fields 
on the same magnitude system with an uncertainty less 
than 0.005 magnitudes (internal error). 

3.3. Artificial star tests 

In order to obtain the LF, we estimated the complete- 
ness of our samples. The completeness corrections have 
been determined by standard artificial-star experiments 
on our data. For each field, we performed five indepen- 
dent experiments adding on the original frames as large a 
number as possible of artificial stars in a spatial grid as 
in Piotto & Zoccali (?). We performed the photometry of 
the frames with the added artificial stars following all the 
steps and using the same parameters as for the original 
images. We converted the new instrumental magnitudes 
to the standard photometric system using the coefficients 
of the calibration equations of Section 2.2. 

In the LFs we include only the values for which the 
completeness corrections, defined as the ratio between the 
number of found artificial stars to that of the original 
added ones, were 50% or higher. This sets the limiting 
magnitude to V = 24. 

4. The Color-Magnitude Diagrams 

We obtained the V and I photometry with the aim of de- 
termining the LF (and MF) of the main sequence stars in 
NGC 4815. We used two colors as we wanted also to de- 
rive the CMD, which allows to: i) discriminate stars from 
false detections; ii) better identify the cluster population, 
in particular the MS stars. 

The CMD for NGC 4815 is shown in Fig. [j] whereas the 
CMDs for the four nearby fields are presented in Fig. 
|2|. Clearly, the CMD of the cluster field does not differ 
much from those of the surrounding Galactic disk. This 
means that NGC 4815 is actually badly contaminated by 
background/foreground stars, as expected from its partic- 
ular position in the Galactic thin disk. The cluster field 
(Fig. 1) is much more populated than the neighboring 



fields (Fig. 2), confirming that NGC 4815 is a real open 
cluster (see also the discussion in Section 5). Our photom- 
etry reaches V « 25 at the base of the MS. 
The MS extends almost vertically for more than 10 mag- 
nitudes, from V = 15 to V=25, although the completeness 
analysis prevents us from using stars dimmer than V= 24. 
The width of the main sequence gets greater at larger mag- 
nitudes. This is partly due to the increasing photometric 
errors at increasing magnitudes (see Fig. 1). However the 
MS is much wider than expected simply from photometric 
errors. Between the various causes which concur to enlarg- 
ing the natural MS width we can envisage the presence of 
unresolved binary stars, common in open clusters, and a 
possible spread in metals. The red vertical sequence, which 
detaches from the MS at about V = 22.5 represents the 
giant branch of the Galactic disk population. 
At odd with the previous photometry, the present data 
(see Fig. 1) do not allow us to recognize clearly the pres- 
ence of a red giant clump, better visible in the BV pho- 
tometry of Carraro and Ortolani (?) at V ~ 14.3. This 
is due to the saturation of many bright stars above V = 
14.5, which were not of interest for the purposes of this 
work. 

5. Star counts and cluster size 

The aim of this section is to obtain the surface density 
distribution of NGC 4815, and derive the cluster size and 
extension in the magnitude space by means of star counts. 
The cluster radius is indeed one of the most important 
cluster parameters, useful (together with cluster mass) for 
the determination of cluster dynamical parameters. Star 
counts allow to determine statistical properties of clusters 
(as visible star condensations) with respect to the sur- 
rounding stellar background. 

Star counts have been performed following the meth- 
ods described in (?), (?) and (?) by evaluating the func- 
tions N(r) and N(V) (N(I)). N(r) represents star number 
in the circle of radius r and N(V) is star number up to 
apparent magnitude V, i.e. integral luminosity function. 
In order to obtain N(r) we divided the cluster field (and 
comparing fields) in concentric rings 40 pixels (11 arcsec) 
wide and counted the number of stars falling inside each 
ring. Then N(r) is obtained by summing the correspond- 
ing ring star numbers. Similarly, in order to obtain N(V), 
we divided stars in the various fields in magnitude bins 
0.2™ wide and summed up corresponding bin star num- 
bers. Star counts with N(r) give us cluster radius and 
cluster star number, star counts with N(V) give us the 
magnitude of faintest cluster stars (up to the limiting mag- 
nitude). 

Due to the saturation of the brightest stars (V < 14.5), 
some regions of the cluster cannot be used to perform star 
counts. In these cases, for N(r) determination we used the 
stars located in the same regions from the work of Carraro 
& Ortolani (1994). 
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Fig. 2. CMDs of the observed nearby fields. Each field is centered at 30' from the center of NGC 4815. The southern 
field do not have deep exposures, and has not been used in the paper. 



A smoothed surface density F(r) is then obtained by 
differentiating numerically the best fit least squared poly- 
nomial approximation for N(r) according to the equation: 



N(r) = 2tt / r'F{r')dr' 
Jo 



(1) 



(sec Sclcznev (1994), for details). The results are shown in 
Fig. 3, where the logarithm of the surface density profile is 
plotted as a function of the limiting magnitude. From this 
figure it is evident that the cluster significantly emerges 
from the background only for V less than 21. Below this 
value NGC 4815 is completely confused with the Galactic 
field. Moreover we note a significant density decrease at 
the cluster position. 



Integral distribution functions are then obtained both 
for cluster field and for comparing fields. Let us consider 
below a theoretical integral distribution function JVj(x), 
where x can be r, V and /, and i — indicates cluster 
field, while i = 1, 2,3, ... indicate comparing fields. 

Following (?) and (?) we consider the normalized de- 
viation z of Ni (x) from the corresponding average number 
of field stars: 



Zi(x) 



N,(x) - N(x) 



a[x) 



(i = 0, !,...,*;), 



(2) 
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where 



i=l 

and 



a 2 (x) = 



k-1 



k 

E 

i=l 



[Ni(x) - N(x) 



(3) 



(4) 



The normalized deviations Zi are connected with the prob- 
ability that the i-th field significantly deviates from other 
fields. Then z values are usually called "significance lev- 
els". 

There are two criteria to determine cluster boundaries 
in x-space. The first criterion is as follows. One has to look 
for a value x c such that (Nq(x) — N(x)) — N c = const 
for all x > x c and zq(x) > Zi(x) (i — 1,2,3,...) for all 
x < x c . In this case we can conclude that the cluster 
has the boundary at the point x — x c and N c is cluster 
star number, and we can estimate an error for x c and iV c 
(Danilov et al. 1986). 

The second criterion is as follows. The cluster has no 
boundary in the sense of the first criterion, but for x > x c i 
zq(x) < Zi(x) for some i, i.e. i-th field deviates from mean 
N(x) more significantly than cluster field. This situation 
means that we can determine only a lower estimate of the 
cluster boundary, i.e. x c > x C 2 and N c > N{x C 2). The 
reason is related to large-scale fluctuations of integral dis- 
tribution function in the comparing fields. Unfortunately 
this situation frequently occurs for open star clusters lo- 
cated in a strongly fluctuating background. 

The star counts in V and / spaces have been per- 
formed both for central part of NGC 4815 (circle with 
r ~ 2 arcmin in the central field, see Table 2) with twelve 
comparing fields (four circles in each East, West and North 
fields) and for all cluster field (i.e. the entire central field) 
with three comparing fields (East, West and North). 

The star counts for the central part of the cluster 
confirm that it is less populated by faint stars (V > 
lgm _ 20"') than the other fields (see Fig.3). We obtained 
cluster boundary V c = 20.8 m ±0.2 m following the first cri- 
terion (see above). The trend of zq(V) is shown in Fig. 4, 
where it is compared with mean z for North, West and 
East fields. The luminosity function of cluster stars (see 
following section), that is differential distribution func- 
tion, has negative value for V > V c . Presumably, the rea- 
son is the higher light absorption along line of sight beyond 
the central part of the cluster. 

By performing star counts in V and / for all the cen- 
tral field we did not obtain the cluster boundary accord- 
ing to first criterion. In the / passband No {I) significantly 
deviates from mean N(I) up to I « 22 m , but in the V 
No(V) significantly deviates from mean N(V) only up to 
V w 20 m . However (N (V) - N(V)) and (JV (J) - N{I)) 
increase for much fainter magnitudes, and the correspond- 
ing luminosity functions are positive (see Fig. 5 and 6 in 
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Fig. 3. Radial surface density profiles of NGC 4815 for 6 
different limiting magnitudes. Surface density is in stars 
per arcsec 2 . 



the following section). This contradictory result is due to 
the very complicate character of the stellar background in 
the vicinity of NGC 4815 and beyond it. It means that 
mean integral distribution function N(V) in comparing 
fields increases more slowly than iVo(V), but the disper- 
sion a(V) is too large. 

Because the cluster center is slightly shifted with re- 
spect to the center of the central field, star counts in r- 
space can be carried out in the entire circle only to radius 
r = 880 pixel (4 arcmin). Due to this reason star counts 
were made also in sector of 199° to radius r = 1080 pixel 
(4.9 arcmin). The results of star counts for different lim- 
iting magnitudes are listed in Table 4. 
First column contains limiting magnitude, second and 
third columns contain cluster radius and cluster star num- 
ber from counts in the entire circle, fourth and fifth 
columns contain the same values from counts in 199° sec- 
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Entire circle {r max = 


880 pixel) 


Sector 199° (r max = 1080 pixel) 




R c 


N c 


R c N c 


16 m 


800 ± 71 


89 ±6 


800 ±55 89 ± 2 


17 m 


800 ± 79 


219 ±13 


800 ± 40 261 ± 3 


18 m 


800 ± 56 


361 ± 11 


840 ± 56 454 ± 5 


igm 


800 ± 64 


447 ± 22 


> 1080 > 637 


20 m 


> 680 


> 319 


> 1080 > 599 




Cluster does not deviate from comparing 


fields significantly, N c > 


22 m 


Cluster does not deviate from comparing 


fields significantly, N c < 


23 m 


Cluster does not deviate from comparing 


fields significantly, N c < 


24 m 


Cluster significantly deviates from comparing fields as minimum, N c < 


24.7 m 


Cluster significantly deviates from comparing fields as minimum, N c < 0. For r 6 
[800, 880] N c > and for r 6 [840, 880] cluster significantly deviates from comparing 




fields as maximum. 
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Fig. 4. Significance level for cluster and surrounding 
fields. See text for details. 

tor. Star number in fifth column is obtained as result of 
star counts multiplied by (360/199). 

The comparison of results of star counts in the entire 
circle and in the sector shows that the asymmetry of the 
cluster increases at increasing limiting magnitude. Begin- 
ning from Vum = 21 m the cluster disappears as a sta- 
tistically significant stellar density fluctuation, and below 
Vum = 24 m the cluster region appears as stellar density 
minimum - statistically significant comparing with three 
fields. The most probable reason is the presence of absorb- 
ing cloud just beyond the cluster. 



Star counts with stars up to Vu m = 19 m gave us the 
estimate of cluster radius, which turns out to be 3.6 ± 0.3 
arcmin, significantly lower than 4.6 arcmin estimated by 
Carraro & Ortolani (1994). However, result of Carraro & 
Ortolani (1994) can be regarded only as preliminary, be- 
cause cluster radius was determined only by surface den- 
sity profile as distance to a point where the slope of this 
curve changes, being the slope negative along all the pro- 
file. It is difficult to compare our results with previous 
results directly, because now we are using a quantitative 
statistical method. Unfortunately, the cluster radius of 
3.6 arcmin is very close to the limit of investigated field, 
and we can not regard this result as definitive. Further- 
more, 3.6 arcmin corresponds to a radius of 2.6 pc at the 
distance of 2.5 kpc. This value lies in the lower limit of 
cluster radii interval according to (?). Then it is quite 
possible that we have determined the radius of the cluster 
core. 

It is necessary to note that the central field contains 
much more stars than any comparing field (sec Fig.l and 
Fig. 2). This excess in star number is mainly due to the pe- 
riphery of the central field which is very rich in faint stars. 
Unfortunately our star counts can not unambiguously an- 
swer the question - is the cluster really disappearing at 
Vum = 21™ and has radius of about 3.6 arcmin or faint 
stars in the outskirts of central field belong to cluster and 
form its halo ? 

Firstly, star counts in wider fields are necessary (both 
for cluster and comparing fields). Desirable radius of in- 
vestigated field is of about 10-15 pc (i.e. 14-21 arcmin). 
Secondly, it is necessary to study color extinction for stars 
of different magnitudes, i.e. to observe in three colors. 
On the basis of the previous discussion, we are going to 
estimate the luminosity function for all the central field in 
the magnitude range V e [14 m ; 24™]. 
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6. Luminosity Function 

The luminosity function of NGC 4815 was obtained as 
follows. In the CMDs, two lines, to the left and right 
of the main sequence, defined by the two relations V = 
6.8(V -I) + 12 and V = 6.8(V - I) + 3.84, were used in 
order to select only the main sequence stars (Fig. 1). For 
each field, an adaptive kernel estimator was then applied 
to obtain a continuous LF. The LF was also derived by 
the standard histogram technique. A detailed description 
of the application of the kernel function method for esti- 
mating distribution functions can be found in (?), and in 
(?) and (?). 

The advantage of this method is twofold: 1) first of all 
it is not necessary to define a bin, whose size usually is 
arbitrary and does not allow to detect unambiguously all 
the statistically significant features present in the LF; the 
estimate of the LF by a kernel method with its confidence 
interval is an estimate of an unknown density of a stellar 
population, while the histogram alone gives only the dis- 
tribution of a specific sample; 2) the form of a histogram 

- crucial for the determination of the LF slope - depends 
both on the bin width and on the initial bin location. The 
distribution function provided by the kernel estimate only 
depends on the kernel width. Moreover criteria have been 
developed to fix the kernel width with high degree of con- 
fidence (?). 

The adaptive kernel estimate of the density distribu- 
tion in one-dimensional case is calculated by the formula: 

/w = (5 ) 

2—1 X ' 

where n is the size of the observing sample Xi...X n of the 
random variable x, K is the kernel function which must 
satisfy the condition: 

/jrw*-!, (6) 

— OO 

and Aj are the local window width factors. The adopted 
kernel has the form 

when \x\ < \/5, and zero elsewhere. 

Basically eq.|| spreads each observed value Xi of the 
random variable x over an interval (window) centered on 

X-i with the probability density K ( x ^\ Ci J ■ The kernel 
width h\i can differ from point to point. The advantage 
of having an adaptive kernel is that it is easy to properly 
sample the tails of the distribution, just playing with the 
kernel width, which is set larger in low-density regions 
automatically (see (?)). Parameter h was set to be 0.3 m . 



i — r 
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Fig. 5. The V LFs obtained with both the standard his- 
togram technique (triangles) and with an adaptive kernel 
estimator (circles). The dotted curves are the confidence 
intervals for the adaptive kernel estimator. Note the large 
error bars below V — 21, when the clusters starts to vanish 
into the Galactic field. 

To determine the confidence interval for the resulting 
LF, we used the method described in Seleznev et al. (2000) 
and references therein. 

For each field, the LFs have been corrected for incom- 
pleteness. The LFs of the external fields were used to sub- 
tract the field-star contamination from the LF of the cen- 
tral field. 

To this end we defined fi(V) as the density of the dis- 
tribution of the stellar V magnitudes, where i is the num- 
ber of the field (here we assume that i = corresponds 
to the cluster field or central field and i = 1,2,3 corre- 
spond to the external fields). The luminosity function of 
the cluster stars is given by 

k 

/(to = /o(to-~ £/i(n (8) 

<=1 

The same procedure was repeated to get the I LF. 
The V and I LFs obtained both with the histogram and 
the adaptive kernel techniques for all the central field are 
shown in Figs. 5 and 6. Note the large uncertainties below 
V = 21 (7 ~ 19), where fluctuations of stellar background 
density play a significant role. LFs obtained by different 
methods are close each other due to large volume of the 
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Fig. 6. As for Fig. 5, but for the I filter. 

sample. Note that the stellar content below V = 20 (I = 
18) is due to outskirts of the central field (see previous 
section) . 

7. Mass Function 

The LFs can be transformed into MFs using a mass- 
luminosity relation (MLR). Since we cannot obtain an 
empirical transformation, we must rely on the theoretical 
models. Therefore, we used the ZAMS isochrones by 
Girardi et al. (2000) with a metallicity Z=0.008. The 
distance modulus and the reddening that result from the 
fit with the model are in agreement, within the errors, 
with the values found by Carraro and Ortolani (1994). 
We adopted E(B - V) = 0.70 and (m - M) v = 14.10. 
From these values we get (to — M)j = 13.20 (see Table 
1). Fig. 7 shows the MF from the V LFs obtained by the 
adaptive kernel estimator and the standard histogram 
technique. The two MFs are rather similar. There is also 
a substantial agreement between the MFs from the V and 
I LFs (Fig. 8), showing an internal consistency of both 
the models and the method. 

By looking at Fig. 7, it appears that there is a dip in 
MF at to ~ 0.9M. . This dip in MF is a the consequence 
of the corresponding dips in LFs at V ~ 19. 5 m — 20™ 1 
and / ~ 17.5'" - 18 m (see Fig.s 5 and 6). To estimate 
the statistical signifincance of this dip, we concentrate on 
"Adaptive kernel estimator" method, which is is more 



Fig. 7. Comparison of the MFs of NGC 4815 from V LFs 
derived by the adaptive kernel estimator and the standard 
histogram technique. The dotted curves are the confidence 
intervals for the adaptive kernel estimator. 

sensitive to peculiar properties of LF and MF than the 
histogram technique (see previous sections). 
We can use the ratio r\ of the depth of this minimum 
(relatively to the mean smooth MF level) to the dis- 
persion of MF estimate in the vicinity of MF minimum 
as a measure of the statistical significance of the MF 
minimum. As for the MF minimum at to ~ 0.9M. we 
obtain rj ~ 1.98 from the "adaptive kernel" MF estimate, 
which corresponds to a reality probability of ~ 95%. On 
the other hand, 77 ~ 0.78 comes out of the "histogram" 
MF estimate, which corresponds to reality probability of 
~ 56%. Then we can conclude that this dip of MF seems 
to be statistically significant. 

Finally, it is possible to obtain an estimate of the clus- 
ter mass by integrating the luminosity function. In detail, 

Mv max 

M=<m>xN = N J m(M v )(f)(M v )dM v (9) 

M Vmm 

where My is absolute stellar magnitude My — V — 
(to — M)y, 4>(My) is the normalized luminosity function 

4>(M V ) = f(M v )/N, (10) 

N is the total number of stars, m(My) is the adopted 
mass-luminosity relation. The luminosity function and 
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Fig. 8. Comparison of the MFs of NGC 4815 from the V 
and / LFs derived by the standard histogram technique. 

mass-luminosity relation have been approximated by a 
spline. Limits in the integral are My min = —0.1 and 
My max — 9.9. The integral turns out to be 

M ~ 1620 ± 690M© 

for kernel LF estimate and 

M ~ 1820 ± 650M Q 

for histogram LF estimate. 

This mass estimate refers to total central field. Unfor- 
tunately it is very difficult to evaluate what part of total 
cluster mass is constituting by this estimate. It depends 
on both unknown cluster size and unknown upper limit of 
magnitude. If we assume that NGC 4815 has a limiting 
magnitude of V c = 20. 8 m , we get the mass estimate 

M ~ 880 ± 230M Q 

The error is determined from mass estimate for upper and 
lower boundaries of confidence interval for LF (see Fig. 5). 

8. Discussion and conclusions 

The MF can not be represented by a unique power law. 
In the mass range 2.5 < m/ra© < 0.8, the PDMF is re- 
produced by a power-law with a = 3.1 ± 0.3 (a — 2.35 
for the Salpeter MF in this notation). In the same mass 



Fig. 9. Comparison of the NGC 4815 MF with other open 
cluster MFs. For display purposes the MFs are shifted by 
an arbitrary constant. 

range the slope of the MF from the I LF is a — 2.9 ± 0.3. 
Below O.8m , the MF is much flatter (a = -0.68 ±0.49). 
It is difficult to correctly interpret the PDMF in this low 
mass regime. However we believe that the flattening of the 
PDMF below O.8m0 can be due to basically two effects. 
First of all, we have already noted how the cluster van- 
ishes (but not disappears) into the Galactic field below 
V = 21, i.e. m = O.8m0. This is likely originated by irreg- 
ularities in the stellar background. As a consequence low 
mass stars are preferentially missed due to absorption. 
Secondly, we cannot however exclude cluster dynamical 
evolution, which causes a loss of preferentially low mass 
stars. Also, stripping of low mass stars can occur along 
the cluster orbit, due to the interaction of the cluster with 
the dense Galactic field. NGC 4815 is indeed located close 
to the Sagittarius-Carina spiral arm, rich in HII regions 
and molecular clouds. Estimating the relative importance 
of these two effects is very difficult and goes beyond the 
aim of this paper. 

Therefore, the slope of the MF below 0.8 tuq can not be 
considered representative of the cluster IMF. 

On the other hand, whether the PDMF above 0.8 
is representative of the IMF is a matter which needs a 
careful model of the cluster too. 

A recent review on the IMF derived from the star counts 
is in Scalo (?). From the comparison of several studies, 
he found that the spread in the slope at all masses above 
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about lm is too large for adopting some average value, 
and in particular, there is no evidence for a clustering 
of points around the Salpeter value (2.35). Scalo (?) con- 
cluded that either the systematic uncertainties are so large 
that the IMF cannot yet be estimated, or that there are 
real and significant variations of the IMF slope. For sure, 
the dynamical evolution, which differs from cluster to clus- 
ter, can modify the MF slope, and also this effect con- 
tributes to the MF slope dispersion. It has emerged that 
the MF can not be reproduced with a single power low: at 
low-masses the IMF from the field stars, open and glob- 
ular clusters (?) is relatively flat, and it becomes steeper 
above about lm©. Also NGC 4815 seems to show this 
break in the MF shape at m ~ lm Q (Fig. [?]), though it 
is not possible to disentangle the real differences in the 
IMF from the effects of the dynamical evolution. A final 
point is the effect of unresolved binaries and mass segrega- 
tion. As discussed by Scalo (1998, and references therein), 
unresolved binaries and mass segregation can only flatten 
the apparent IMF derived by not accounting for these, so 
the real IMF must be steeper than the quoted slope by an 
unknown amount. 

In Fig. |, the MF of NGC 4815 is compared with the 
MFs for other 15 open clusters. Most of them cover masses 
above lm©. We have derived the MFs directly from the 
published LFs using the same set of isochrones from Gi- 
rardi et al. (?) used to get the NGC 4815 MF. The MFs in 
Fig. ^ are shifted by arbitrary constants for display pur- 
poses. The MF of NGC 4815 is one of the most extended 
(in mass) MFs obtained so far for an open cluster. Only 
for the better known and nearby open clusters (Pleiades 
and M35), and for objects observed with HST (NGC 2516, 
NGC 2477, and NGC 2420) do the MFs reach ~ 0.1 m©. 

Fig. |l(^ shows the slopes a of the MFs plotted in Fig. || 
as a function of the mean mass. The slopes have been 
obtained by fitting a power law in the m < lm© and m > 
Ira© mass ranges. Clusters with ages greater and less than 
or equal to 0.5 Gyr are plotted with different symbols. As 
noted by Scalo, the spread in slope is very large. However, 
it is interesting to note that all the clusters with ages 
< 0.5 Gyr have very similar slopes (a = 2.25 ± 0.88) for 
m > lm©. 

NGC 4815, which has an age of 0.5 Gyr (Carraro and 
Ortolani (?)), has a MF slope fully consistent with this 
value. The MFs of the clusters older than 0.5 Gyr have 
a similar or flatter slope. This may indicate that their 
PDMF may have not or may have been changed by the 
dynamical evolution which has had more time to operate. 
A dynamical model of each of these clusters is necessary to 
confirm this interpretation of Fig. [Io| At the moment, the 
similarity of the MF slopes for all the youngest clusters in 
Fig. [l^ can be interpreted as an evidence of a similarity of 
the IMF, which seems to not depend much on metallicity, 
stellar density, or galactocentric radius. 

The interpretation of the MF slopes below 0.8m© is 
more difficult, both because there are less (and more un- 



Scalo, 1998 

Salpeter, 1959 
Scalo, 1998 



NGC 1758 



• ages > 0.5 Gyr 
o ages < 0.5 Gyr 

A NGC 481 5 

_j , I , 

-0.5 



0.5 



log m/m E 



Fig. 10. The MF slopes a are plotted as a function of the 
mean fitting mass in the two mass ranges: m > lm© and 
777 < lm Q . Filled circles represent clusters older than 0.5 
Gyr, whereas open cirles stand for younger clusters. Note 
the similarity of the MF slopes for the youngest clusters. 

certain) data, and because we expect that in this mass 
regime the dynamical evolution should be more effective 
in altering the IMF. Still, all the MFs seems to have the 
same flat general trend. We also note that these slopes 
represent an upper limit for the slopes of the MF of the 
Galactic globular clusters (Piotto and Zoccali (?)). 
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